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Abstract 


Adaptive grid generation is an active research topic for numerical so¬ 
lution of differential equations. In this paper, we propose a variational 
method which generates transformations with prescribed Jacobian deter¬ 
minant and curl. Then we use this transformation to achieve adaptive 
grid generation task. 
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1 Introduction 

The deformation method 14],[5],[7] has its origin in differential geometry [I]. 
Its main theoretical advantages is that the grid transformation generated has 
prescribed Jacobian determinant. Its implementation is based on calculating 
a velocity field from a differential system consisting of the divergence operator 
and the curl operator. In the special case when the curl vector is set to zero, the 
differential system reduces to three (decoupled) scalar Poisson equations in 3D. 
Among existing methods are several interesting variational methods 
which enjoy clear physical interpretations and are straight forward to implement. 
The problem we study in this paper is to construct numerically a transformation 
T on a two dimensional domain D, such that, the Jacobian determinant of 
T, namely J(T), equals a prescribed monitor function /o, and the curl of T, 
citrZ(T), equals a prescribed monitor function go: 


J(T(x)) = /o(x) 


( 1 ) 


curl{T{x)) = go{x) 


We define the — norm of the difference between J(T) and /o, curl{T) and 
go as 
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where a > 0 is a weight parameter. 

Our goal is to construct the desired transformation T by minimizing ssd, with 
respect to some control function f under the following constraints: 

T(x) = X + u(x) (3) 

Where the displacement u satisfies the Poisson’s equation 

Au = / (4) 

with fixed boundary conditions. Given arbitrary 5f, we have Jf = JAu = A((5u). 
So Ju is governed by (5f through a Poisson’s equation. 

In order to use the gradient descent method to minimize ssd with respect to 
control function f, we shall derive the gradient based on the fixed Dirichlet 
boundary condition. Then we will demonstrate the method by examples of 
recovering some known transformations. 


2 Derivation of the Gradient 

We focus on the case where the grid nodes on the boundary are fixed. Since 
the boundary nodes are fixed, the method enjoys convenience in computation 
and elegant derivation in mathematics. It turns ont the general cases where 
boundary nodes are moved, can be reduced to this fixed boundary case, see 
Remark I at the end of this section. 

Now we let T = (Ti,T 2 ),u = (■ui,'U 2 ),f = (/i,/ 2 ), and take a = 1 for the 
simplicity of presentation. For the general cases, see Remark 2 at the end of 
this section. So, 

ssd = i JJjiJiT) - fo f + (curliT) - go f]dA 

and the displacement u satisfies now, 

( Aui = /i 

< Au 2 = /2 (5) 

y ui,U 2 =0 on dfl 

By variational calculus, we have 

Sssd = //fj(J(T) - fo)SJ{T) + icurl{T) - go)Scurl{T)dA 

= //f 2 [(./(T) - fo) 6 iTi,T 2 y - TiyT 2 oo) + (curKT) - go)SiT 2 ^ - Tiy)]dA 

= JQiJiT) - fo){STi,T2y + TiJT2y - 5TiyT2^ - Tiy6T2^) 

+ {curl{T) - go){5T2x - 5Tiy)]dA 

From (3), we know (5T = (5u, or 5Ti = Sui,ST 2 = Su 2 , so STix = Suix,STiy = 
Suiy, z = 1,2, plug back in the above equation, we can continue as: 

Sssd = jj [{J{T) - fo){ 5 uixT 2 y + Tix 5 u 2 y- 5 uiyT 2 x-Tiy 6 u 2 x) 

+ {curl(T) - go){5u2x - 5uiy)]dA 
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Let P = (J^(T) — /o), Q = {curl{T) — go), we have 


Sssd = 



[P{6ui,^T2y + TiooSu2y 


- duiyT2x - TiySu2x) + QiSu2x - Suiy)]dA 


Now, notice that VSui = {5uix,5uiy),i = 1,2 


5ssd = 



[P{{T2y, -T2 x) • V(5mi + {-Tiy, Tlx) ■ V(5u2) 


+g((0, -1) • W5ui + (1, 0) • V5u2)]dA 


Let 

r -ai =P{T2y,-T2x) + Q{0,-l) 

\ -a2 =P{-Tiy,Tix)+Q{l,0) 

(6) 

we get 

Sssd = (—ai • ySui — a 2 ■ VSu 2 )dA 

JJQ 



Before further proof, we need to list 2 corollaries of divergence theorem which 
will be used later. 

Corollary 1. By Applying the divergence theorem to the product of a scalar 
function v and a vector field u on a 2 D domain ft, we have 

// {u-Vv+ v{V ■ u))dA= / {vu-n)dS 
JjQ, Jon 

where n is the normal vector. 

Corollary 2. By Applying the divergence theorem to the product of a scalar 
function v and a vector field Vw on a 2 D domain fl, we have 

// {Vu ■ Vv + v{Au))dA = / {vVu-n)dS 
JJQ. JdQ 

where n is the normal vector. 

Notice in Poisson equation ( 5 ), we have Ui = 0 on dD,, so 5 ui = 0 on dD.. By 
Corollary 1 we can continue as 

Sssd = ff^(—ai ■ VSui — a2 • V6u2)dA 

= //q((V • ai)( 5 Mi + (V • a2)6u2)dA - fg^(Sui(ai ■ n) + 6u2{a2 ■ n))dS 
= ■ ai)6ui + (V • a2)6u2)dA 

Now Let g = (gi, (72) satisfies 


A^i = V • ai 
Ag2 = V • a2 

=0 on dVt 


( 7 ) 
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then we have now, 


Sssd = / / (Agi5i6i + Ag 2 Su 2 )dA 
J Jn 


By the Corollary 2, for t = 1, 2 


/ JjQ{'^ 9 i-'^Sui + SuiAg^)dA = Jg^{6uiVgi ■ n)dS =0 
I + giA6ui)dA = \^{giV6ui ■ n)dS = 0 


along with (5), finally we will have 


5ssd = ff^{AgiSui + Ag 2 Su 2 )dA 

= //n(5i^(^Mi) + g2S{Au2))dA 

= //n (51^/1 + 92Sf2)dA 


which is 


dssd 


1,2 


( 8 ) 


Remark 1. At the end of proof, we want to point out that if we change the 
identity map x to any given transformation T*(x) in (3), the derivation above 
is still right. Generally speaking, (3) can he replaced by 


T(x) = T*(x) + ^x) (3*) 

where T* (x) is any given transformation. 

Remark 2. For arbitrary given weight parameter a in (2), the derivation above 
is still right by letting Q = a{curl{T) — go) 


3 Implementation and Experiments 

3.1 Implementation by the Gradient Descent method 

Now, we can implement our method by gradient descent method as follows: 

1. Initialize T = id, u = 0, / = 0 

2. Compute ai,a 2 by (6), then solve Poisson’s equation (7) to get ( 71,(72 

3. Update / by fi^new = ft,old — 9 i ^ tstep, where tstep is an optimization 
parameter 

4. Solve Poisson’s equation (5) to get ui,U 2 

5. Update T by (3) 

6. Back to 2, iterate until a preset tolerance or a preset number of iteration 
steps is reached 
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3.2 Experiments of Recovering Transformations 

For the purpose of testing the accuracy of our method, and discovering more 
details inside, we design a recovering experiment. By the derivation, if given a 
known transformation Tq, we can compute its Jacobian determinant and curl 
as our fo,go, namely foix) = J{To{x)), go{x) = curl{To{x)), then the desired 
transformation T constructed here should equal Tq. The following figures and 
tables illustrate corresponding result. In order to compare, we also give the cases 
when only consider the Jacobian term in our method. We can see clearly the 
curl plays an important role in grid transformation, and with both the Jacobian 
and curl term here in our algorithm, we can recover a transformation, which 
proves the algorithm is reliable and accurate in grid generation. 

3.2.1 Recover a transformation Tq with fixed boundary nodes 

Given a nonlinear transformation Tq from the square [1,65] x [1,65] to itself, 
we want to reconstruct Tq from its Jacobian determinant and curl. We assume 
that the boundary nodes of Tq coincide with the nodes of the background grid. 
We have the following results for a = 1,0.1,10 respectively. And for different 
values of a, there is no significant difference. 
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(a) Only Jacobian 



(b) Only Jacobian, enlarged (c) Only Jacobian, enlarged 
view of left rectangle view of right rectangle 



(d) Jacobian and Curl 



(e) Jacobian and Curl, en- (f) Jacobian and Curl, en¬ 
larged view of left rectangle larged view of right rectangle 


Figure 1: Experiment 1.1, a = 1, 65 x 65 grid nodes, 2000 iteration steps, 
9.12sec by Matlab on a Dell 780 PC. The black star dots * represent To, and red dots 
■ represent constructed T 



Only Jacobian 

Jacobian and Curl 

ssdj^ 

0.2640 

0.3301 

ssd 

4.8134 

0.6377 

maximal distance^ 

0.3444 

0.0706 

average distance 

0.1109 

0.0177 

maximal angle difference^ 

25.8343 

15.7230 

average angle difference 

2.6132 

1.7481 


^ ssdj means only consider Jacobian part in ssd, ^distance means 
||T(x) — To(x)||2,x € n, ^angle difference means the corresponding 
angle differences in every corresponding quadrilateral 


Table 1: Comparison of Experiment 1.1 
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(a) Only Jacobian 



(b) Only Jacobian, enlarged (c) Only Jacobian, enlarged 
view of left rectangle view of right rectangle 



(d) Jacobian and Curl 



(e) Jacobian and Curl, en- (f) Jacobian and Curl, en¬ 
larged view of left rectangle larged view of right rectangle 


Figure 2: Experiment 1.2, a = 0.1, 65 x 65 grid nodes, 2000 iteration steps, The 

black star dots * represent Tq, and red dots • represent constructed T 



Only Jacobian 

Jacobian and Curl 

ssdj 

0.2640 

0.4783 

ssd 

4.8134 

1.6869 

maximal distance 

0.3444 

0.0757 

average distance 

0.1109 

0.0198 

maximal angle difference 

25.8343 

16.5480 

average angle difference 

2.6132 

1.8057 


Table 2: Comparison of Experiment 1.2 
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(a) Only Jacobian 



(b) Only Jacobian, enlarged (c) Only Jacobian, enlarged 
view of left rectangle view of right rectangle 



(d) Jacobian and Curl 



(e) Jacobian and Curl, en- (f) Jacobian and Curl, en¬ 
larged view of left rectangle larged view of right rectangle 


Figure 3: Experiment 1.3, a = 10, 65 x 65 grid nodes, 2000 iteration steps. The 

black star dots * represent Tq, and red dots ■ represent constructed T 



Only Jacobian 

Jacobian and Curl 

ssdj 

0.2640 

1.4442 

ssd 

4.8134 

1.8762 

maximal distance 

0.3444 

0.0797 

average distance 

0.1109 

0.0203 

maximal angle difference 

25.8343 

16.6378 

average angle difference 

2.6132 

1.9603 


Table 3: Comparison of Experiment 1.3 


3.2.2 Recover a transformation Tq with moving boundary nodes 

Now we assume that the transformation Tq has different boundary nodes com¬ 
pared to the background grid. We can use the Remark 1 above to find a bound¬ 
ary match transformation T* first, then turn back to the fixed boundary nodes 
case. The following figure 4 shows a experiment for moving boundary nodes 


case. 











































































































































(d) Jacobian and Curl 




(e) Jacobian and Curl, en- (f) Jacobian and Curl, en¬ 
larged view of left rectangle larged view of right rectangle 


Figure 4: Experiment 2, a = 1, 65 x 65 grid nodes, 2000 iteration steps, The 

black star dots * represent Tq, and red dots • represent constructed T 



Only Jacobian 

Jacobian and Curl 

ssdj 

1.7588 

0.4327 

ssd 

107.8698 

1.5451 

maximal distance 

1.0158 

0.1477 

average distance 

0.3316 

0.0204 

maximal angle difference 

43.4929 

16.9158 

average angle difference 

4.7968 

1.2934 


Table 4: Comparison of Experiment 2 


4 Conclusion 

We described a new variational method of grid generation with prescribed Jaco¬ 
bian determinant and curl, the method has very solid and concise mathematical 
foundations, related numerical experiment and results also show the method is 
reliable and accurate. The boundary cases we discussed here almost cover every 
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possible applied case. Also the method can be extended to 3D domain, which 
will be discussed in further papers. 
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